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Abstract. We present a numerical study of the dynamics of the one-dimensional 
Ising model by applying the large-deviation method to describe ensembles of 
dynamical trajectories. In this approach trajectories are classified according to a 
dynamical order parameter and the structure of ensembles of trajectories can be 
understood from the properties of large-deviation functions, which play the role of 
dynamical free-energies. We consider both Glauber and Kawasaki dynamics, and 
also the presence of a magnetic field. For Glauber dynamics in the absence of a 
field we confirm the analytic predictions of Jack and Sollich about the existence of 
critical dynamical, or space-time, phase transitions at critical values of the "counting" 
field s. In the presence of a magnetic field the dynamical phase diagram also 
displays first order transition surfaces. We discuss how these non-equilibrium 
transitions in the Id Ising model relate to the equilibrium ones of the 2c? Ising model. 
For Kawasaki dynamic we find a much simple dynamical phase structure, with 
transitions reminiscent of those seen in kinetically constrained models. 
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1. Introduction 

In this paper we perform a detailed numerical study of the trajectories of the 
one-dimensional Ising model under both Glauber and Kawasaki dynamics. We 
do so by means of the "thermodynamics of trajectories" method introduced in 
references [HIZHHHIl, which in turn is based on Ruelle's thermodynamic formalism for 
dynamical systems (SH6J. The basic idea is to understand dynamical properties of a 
many-body system by treating ensembles of trajectories in a way which is analogous 
to how one understands static or structural properties via the ensemble method 
for configurations (or microstates) in standard equilibrium statistical mechanics. 
In order to classify trajectories one defines a suitable dynamical order parameter, 
one that captures the essence of the physics one is trying to uncover. As we 
will show below, a convenient one is given by the "dynamical activity" |2l ® 0, 
which for a lattice problem like the ones we study in this paper corresponds to the 
number of configuration changes (spin-flips in our case) in a trajectory. Alternatively, 
it could be the time integral of a standard observable, such us the time-integral 
magnetization. One can then study ensembles of trajectories classified by dynamical 
activity, or more conveniently as we show below, by its conjugate "counting" field 
s. For trajectories observed over long-times a large-deviation principle (HI holds 
such that the distributions of the activity over the ensemble of trajectories, or their 
corresponding generating functions, are determined by time independent "large- 
deviation" functions, which play a role analogous to thermodynamic potentials 
(entropies or free-energies) for the dynamics. This approach has proved successful 
recently in uncovering dynamical or "space-time" phase transitions in a variety 
of classical many-body systems, most notably glass models, see e.g. ||3l |9j [TO]. ITT1] . 
In these systems thermodynamics is either trivial, as in kinetically constrained 
models of glasses [12], or not very revealing, as in more realistic liquid models Ili~3i 
but the dynamical phase structure is extremely rich: i.e. complex and highly 
correlated relaxational dynamics can emerge independently of any form of complex 
thermodynamics, and this is clearly revealed by the large-deviation method for 
trajectories. 

So a natural question is the following: what is the dynamical phase structure 
of systems which do display strong thermodynamic features, and is there any 
relationship between thermodynamic and dynamical, or space-time, phases? Recent 
work has addressed this issue in certain simple models, including mean-field Ising 
and Potts models [14] and certain mean-field spin glass models |[T5l[T6| . Specifically, 
and directly related to the work we present here, there is a recent paper by Jack 
and Sollich |i~7f in which they study analytically the dynamical phases of the Id 
Ising model with Glauber dynamics. By mapping the problem to that of a quantum 
Ising model in a transverse field, they managed to obtain the large-deviation 
function that plays the role of a free-energy for trajectories, and showed that in an 
enlarged T — s (temperature /counting-field) space of parameters there are second- 
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order dynamical transitions between active (and paramagnetic) and low-activity 
(and ferromagnetic) dynamical phases — transitions which are not observed under 
standard thermodynamic conditions — with critical properties of the classical 2d Ising 
universality class. 

The current paper builds on Jack-Sollich's results. By means of path sampling 
methods 11181 we generate ensembles of trajectories biases by the counting field s, 
the so-called s-ensembles (9). We confirm Jack-Sollich's predictions for Glauber 
dynamics in the absence of the magnetic field, providing also detailed analysis of 
the scaling properties of the dynamical transitions. We also extend the problem to the 
case of non-zero magnetic field. This allows us to map out in detail the full dynamical 
phase diagram, revealing a host of first- and second-order dynamical transitions, and 
the connection between dynamical and thermodynamical phases. We also consider 
the case of Kawasaki dynamics and show that the dynamical phase structure in this 
case is analogous to that of kinetically constrained models of glasses B3HU. 

This paper is structured as follows: first, in section |2], we will introduce the 
theoretical framework and large deviation formalism, followed by two main sections. 
In section |3] the theoretical and computational results of the Id Ising model using 
Glauber dynamics are discussed with strong focus on the extension of the phase 
diagram through an external magnetic field and the resulting scaling behaviour. 
In the section H] we will look at the dynamic behaviour when applying Kawasaki 
dynamics. Finally our conclusions are given in section [5j 

2. Theory and Methodology 

The objective of this paper is a computational study of the dynamical processes of 
a Id Ising spin chain in the presence and absence of a magnetic field, and within 
Glauber or Kawasaki dynamics. In order to achieve this, we apply what is referred 
to as the s-ensemble in the literature O HI, which will be briefly reviewed in the 
following. In order to study the dynamics, or time evolution of the system of interest, 
an appropriate observable capturing the dynamics needs to be defined, as well as a 
set length of time for which the system is observed. In this way a trajectory or history 
of the system can be defined with a given length in time t Q bs- The trajectory is the 
sequence of configurations, and thus contains information about the instantaneous 
state of the system at each point in time. Any configurational change in the system 
gives rise to a dynamic observable. The most basic one is the number of spin flips in 
a given trajectory. This is referred to as the activity K |]2}|3HZ||. In the thermodynamic 
Gibbs ensemble the probability of observing a specific configuration is governed by 
the partition function Z. For the dynamics an equivalent object can be constructed, 
based on the formalism by Ruelle 10, which gives information about the probability 
of observing a given trajectory. The dynamic partition function is defined as 

Z(s,t ohs ) = (e- sK ) , (1) 
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where K is the activity of the system. We are considering path averages, denoted 
by ( ) , over a set time intervals starting at t — until a later time t = t ohs in an 
equilibrated system. The parameter s is a biasing or counting field conjugate to K 
and hence allows subtle control over the ensemble of trajectories. In the context 
of full counting statistics in electronic transport such parameters are referred as 
"counting fields" [19J. It plays an equivalent role to that of the inverse temperature 
f3 as conjugate field of the energy (or to the magnetic field as a conjugate to the 
magnetization in a magnetic problem) in the Gibbs ensemble. 

In order to give a physical meaning to this space-time analogue of the partition 
function, the large deviation formalism is used. Traditionally the thermodynamic 
limit, as in an infinite system size, is considered. In this case the limit of very large 
observational time t Q b s will be considered. Thus the partition function can be written 
in form of: 

Z(s,t ohs ) ~ e l °^ , (2) 

where the "large-deviation function" El IH |8l plays the role of a free-energy 
function for trajectories. The structure of this function provides information about 
the dynamic phase space behaviour of the system, e.g. it singularities indicate phase 
transitions in ensembles of trajectories. 

In order to analyse the dynamic behaviour we consider a set of time 
extensive observables, well defined within trajectories, such as the time integrated 
magnetization or energy for the case of an Ising magnet. If A is one of these 
observables, its expectation values biased according to the modified path ensemble 
is given by 

(A) - {Ae ' SK)o (3) 

Note that for s = the equilibrium expectation is recovered. Hence a set of time 
intensive variables can be defined, normalized by the choice of observational time 
and system size N. For example for the intensive activity k we have 

*00 = ^ ■ (4) 

This leads to a generalization of the susceptibilities associated to this ensemble, which 
are defined in the usual way 

(K>), - (K)l 

This will further allow the study of the critical behaviour within the dynamics. 

This theoretical framework will be used in order to study a Id Ising spin chain, 
on a lattice of size N, under the presence of an external magnetic field. The general 
Hamiltonian for such a system is given by: 

N N 

H = - J VjVj+i - h Qj, (6) 
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where the spins Oi take values of 1 or —1, depending on their direction, and J > 
is the coupling constant between first neighbours, and takes the value of 1 for our 
purposes. The second sum describes the coupling of individual spins to the external 
magnetic field h. Simulations are done using a Monte Carlo scheme in which periodic 
boundary conditions are considered for the spatial dimension. As usual, the Monte 
Carlo time is taken as one unity when iV attempts of spin flips are performed. The 
trajectory length t obs is predefined. Adjusting the field s, will allow biasing towards 
lower or larger activity and thus samples the tails of the equilibrium distribution. 
Rather than just running different realizations of the equilibrium distribution, we 
use an adaptation of transition path sampling (TPS) with s-ensemble biasing |fT8j, 
similar to the one used to generate s-ensembles in certain glassy models HHH]]. Once 
an equilibrium trajectory is produced, this can be used as blue print for the next 
trajectory. As in TPS only part of the original trajectory is used. A shooting point is 
chosen, from which the system is propagated in time, until the chosen observational 
time is reached. Then, for the original and new trajectory the activity is calculated, 
which is the incremental number of spin flips observed within the length of the 
trajectory. This leads to the acceptance of the newly proposed trajectory according 
to the standard Metropolis criterion: 

P acccpt = min(l,e- sAX ) (7) 

where AK is the difference in activity between the two trajectories. The employed 
TPS method allows the generation of trajectories with similar values of K and 
therefore ensures that new trajectories are likely to be accepted. This method has 
been employed for the generation of data in this paper in an efficient way. 

In terms of the Monte Carlo dynamics used within this work, there were two 
separate computational studies, one using Glauber dynamics and the other Kawasaki 
dynamics. The motivation for the use of Glauber dynamics over the generally 
preferred Metropolis dynamics is theoretically motivated. A recent paper has a 
theoretical solution for the space time behaviour of the Id Ising model with Glauber 
dynamics, employing an ensemble similar to the s-ensemble presented above |T7|| . 
Jack and Sollich refer to this ensemble the ^-ensemble. Rather than biasing towards 
activity an (integrated) energy bias is introduced, where g is the conjugate field to the 
energy. Analogously to the s-ensemble a dynamical free energy <p(g) can be defined 
for this ensemble. Only in terms of the Glauber dynamics a theoretical solution is 
readily available in order to quantify the critical behaviour of the system. Thus, 
the computational model uses a Glauber acceptance criterion, but for the Kawasaki 
dynamics the preferred Metropolis criterion will be used. This will be presented in 
detail in the following two sections. 
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3. Glauber Dynamics 

In this section we will present the space-time phase diagram obtained from our 
computational results using Glauber dynamics. In terms of the evolution of the 
Monte Carlo dynamics, the probability of accepting a spin flip is given according 
to the Glauber acceptance criteria |20] 

-Paccopt = j- — e/3AE (8) 

where (3 is the inverse temperature (3 = l/ksT and the Boltzmann constant &b is 
assumed to be 1. 



3.1. The phase diagram 

Through the temporal evolution a dynamic phase diagram can be constructed in 
temperature and biasing field s. Reference [17] considers a modified path ensemble 
equivalent to the discussed s-ensemble. In the g-ensemble, similar to the dynamic 
free energy ip(s) introduced in section |2l the dynamical free energy 0(g) of the 
biased ensemble of trajectories can be obtained using the large deviations of the 
thermodynamic energy with zero magnetic field. It was further shown, that the 
second derivative of <p(g) diverges to a line of critical points, which can be shown 
to belong to the universality class of the 2d Ising model. An equivalence between 
the activity constrained ensemble (with biased field s) and the energy constrained 
ensemble (with biased field g) can be established, since for the Ising problem we 
are considering biasing the energy is akin to biasing in terms of the escape rate 
of the corresponding continuous-time Markov chain, see e.g. flU. In this model, 
due this equivalence the free energy of the s-ensemble if)(s) shares the same critical 
properties of (f)(g), because both functions only differ in a regular factor. Therefore 
the dependence on s and T of the second order critical line is given by the following 
relationship fl7| 

s(T) = - In [tanh (2J/k B T)} . (9) 

Figured] shows, a schematic representation of the dynamic phase diagram for the 
Ising chain with Glauber dynamics. This phase diagram is defined in a 3d parameter 
space depending on the set of variables {s, T, h}. In particular the line in figure []Ja) 
represents the curve given by Equation dD, which is a line of critical points in the 
plane of h = 0, dividing this plane into two phases: the paramagnetic phase (for 
s < s c ) and the ferromagnetic phase (for s > s c ) |fT7|. 

Figure []Ja) also shows a line for zero magnetic field given by isothermal 
conditions T = T (dashed line). As indicated, this line contains a critical 
point corresponding to s c = s(T ) given by Equation ((9]). Figure [Ob) shows the 
magnetization of the Ising chain as a function of the biasing field s. Here the critical 
point s = s c separates the ferromagnetic phase for s > s c with non zero spontaneous 
magnetization m ^ and, the paramagnetic phase for s < s c with m — 0. Tuning 
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Figure 1. Schematic representation of the phase diagram of the Glauber Id Ising 
model in the s-ensemble with magnetic field h. (a) The phase space is defined 
through three variables T, s, and h. Here the continuous line, into the plane h = 
0, represents the critical points given by the theoretical solution in reference |[T7| 
separating paramagnetic (and active) and ferromagnetic (and low-activity) dynamical 
phases. The dashed line correspond to the conditions T = T = const and h = and 
contains a critical point (CP) marked by a circle, (b) Spontaneous magnetization as 
a function of s at temperature T = Tq and h = 0, (dashed line drawn in (a)). The 
symmetry is broken by the biasing field s for s > s c just at the critical point s = s c . (c) 
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the biasing field s, once the critical point is reached, the symmetry of the system 
is spontaneously broken. The inclusion of the magnetic field in the isothermal 
behaviour gives rise to a third dimension and therefore the plane T = T can be 
considered as shown in the figure Hfc). Here for s < s c the magnetic field acts on a 
paramagnet, while for s > s c the magnetic field acts on a ferromagnet. Figure [Ob) 
and figure QJc) are analogous to those of the 2d Ising model by assuming s — > (3. 
In fact by taking s = s c and using the temperature T (instead of s) as the control 
parameter, the critical behaviour is in complete analogy with the 2d Ising model. In 
the presence of an external magnetic field h we expect the system to exhibit a first 
order phase transition along the dashed line showed in the figure [He) where s > s c . 
This obviously gives rise to an entire surface of critical points for h — below the 
curve of s(T) as showed in figure H{a). In the following, the second order critical 
behaviour as predicted by the theory of Equation © will be studied through Monte 
Carlo simulations as well as the validity of the existence of the first order surface by 
means of the presence of an external magnetic field. 

3.2. Dynamical continuous phase transition 

In this section we study by means of Monte Carlo simulation the critical behaviour 
of the critical points given by the Equation ((9]). As was established in the discussion 
of the section 13.11 describing the figure [H we can use a s-fixed curve or a isothermal 
in order to reach the critical point. Trajectories are generated as described in section |2] 
using a fixed temperature T and with zero magnetic field. From each accepted 
trajectory, of a given value of the biasing parameter s, the set of time extensive 
variables is obtained. These are the activity K as the number of the changes 
of configuration, the integrated energy E = J <obs dt'u(f), and the integrated 
magnetization M = f* obs dt'm(t'), where m(t) = J2? =1 s<(*). 

In order to study critical properties, we need to apply the finite size scaling 
theory. This is a commonly used technique in computational physics plf. As we 
are investigating the scaling behaviour of dynamic properties, the temporal size is 
given by the observational time t D bs- In fact, we can see from the definition given by 
the Equation (|2j) that it is analogous to the spatial system size when the equilibrium 
thermodynamic free energy is defined. Of course, in our simulations we have to use a 
system with finite spatial size, i.e. the number of spins N. We have observed that the 
behaviour of the system is very sensitive with respect to the length of the trajectories 
t bs (much more than to the length N of the chain). In order to have a consistent 
analysis of the finite size effects, we have chosen to study the behaviour considering 
t bs as the scaling parameter, so that the lattice size N is set to fixed value. 

It is worth mentioning that the dynamic evolution of the this system can 
be described with a time evolution operator, which is equivalent to a quantum 
Hamiltonian. From this point of view the length (i bs) of the trajectories is equivalent 
to the inverse of the temperature T^' , as in the associated quantum problem 
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Figure 2. (a) Susceptibility of the energy as a function of T at constant s = 0.036 635 37. 
Vertical dashed line corresponds to T c = 1.00 given by Equation (9). (b) Susceptibility 
at T c — 1.00 versus i b s - The divergence is logarithmic, as in the case of the specific 
heat in the 2d Ising model, (c) Effective critical s, taken as the s at the peak of figure|2] 
(a), against t obs . The extrapolated value, for i obs -t oo, is T c — 1.02 ± 0.02. 



discussed in fl7| (so that t Q bs <^ ). It is well known that there is a mapping 

between the classical 2d Ising model and the quantum Ising chain [22 j. From this 
mapping we know that the inverse of T is equivalent to the size in one spatial 
direction in the 2d Ising model (that is 1/T <-> L). Therefore we conclude that 
tobs in our problem is equivalent to a linear size L in the classical 2d Ising model (i.e. 
t bs L). According to this we assume that the finite size effects are given by the 
same well known relationships used for the thermodynamic 2d Ising model, where 
we now replace L with t Q bs- 

To determinate the critical temperature of a continuous phase transition we can 
use the peak of the susceptibilities. According to the previous discussion we expect 
that susceptibilities have the scaling forms given by ||2T| 

x(Us) = t£f ((T - T c )tlt) , (10) 

where a is the exponent that characterizes the divergence of the susceptibility x an d 
v the length correlation exponent. Furthermore, the location of the peaks define an 
effective size-dependent transition temperature which is expected to vary as 

T c (tobs) = T c (oo) + At2t ■ (11) 

As it is showed in the figure [0 we can use an s-fixed curve or an isothermal one 
in order to reach the critical point, and therefore the scaling given by the Equation 



(|TT|) is expected in terms of T or in terms of s, respectively. 
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Figure |2] shows the simulation results for a fixed value of s and varying the 
temperature T. We have chosen the value of s as s c = 0.036 635 37, which correspond 
to s c at s(T = 1) as given by Equation |9|). We define the generalized susceptibility of 
the energy x e by means of 

= m, - {sr. 

In figure |2{a) we plot x e as a function of the temperatures and for different trajectory 
lengths t obs - We used a lattice size N = 64, and the averages were taken over 5 x 10 6 
realizations of single trajectories. Peaks close the critical point are clearly observed, 
which sharpen, as the system size is increased (i.e. t Q bs)- We have used the location of 
these peaks as the definition of a size-effective "critical temperature" T c (t ohs ). 

Figure |2][b) shows the effective critical temperature as a function of the inverse of 
the trajectory length £ obs . We can see that the behaviour is compatible with Equation 
((TTj) taking an exponent v = 1 corresponding to the 2d Ising universality class. The 
continuous line is given by the linear fit of the data. From this our extrapolated 
estimated value for the critical temperature is T c = 1.02 ± 0.02, and it is consistent 
with the theoretical value given by Equation (0). Also, figure |2tc) shows the value 
of the susceptibility measured in figure |2](a) for T = T c = 1 and different sizes in 
log-linear scale. 

As the dynamical free energy 0(g) must have the same critical properties as the 
thermodynamic free energy f(T) of the 2d Ising model, its second derivative must 
be analogous to the behaviour of the specific heat. In fact it can be shown that 
the second derivative of 4>{g) has a logarithmic divergence when the critical point 
is approximated fT7|. Figure |2{c) shows that a logarithmic behaviour is observed, 
which means that Equation ((T0|) holds with an exponent a = 0, according to the 
previous sentences. 

As we have already indicated, the critical behaviour can also be accessed by 
choosing a constant temperature T and considering different values of the biasing 
field s. In figure |3ta) we plot the susceptibility of the activity %k' given by 
the Equation ©, versus s with fixed T = 1. The observed behaviour is analogous 
to the energy fluctuations. Again the value of the critical s c can be extrapolated 
making use of Equation (|TT|), but in this case using the biasing field s instead of T. 
In figure |3](b) we plot the effective critical s and the linear fit gives us an estimate of 
the critical field, that is s c = 0.036 ± 0.01. Once again this result is consistent with 
v = 1. Also figure |3](c) shows the susceptibility Xk for s = 0.036 635 37 in a log-lineal 
scale. We observe again that the Equation (fTO]) is valid with a logarithmic divergence 
of the susceptibility of the activity, that means a = 0. 

So far a complete equivalence, in the critical behaviour, between the biasing 
field s and the inverse of the temperature 1/T is observed. Also, similar results are 
obtained when studying finite size effects of the fluctuations of K (or E) with respect 
to T (or s). 




Figure 3. (a) Susceptibility of the activity as a function of s at constant temperature 
T = 1. Vertical dashed line corresponds to s = s c = 0.036 635 37 given by Equation 
©. (b) Convergence of the effective critical s, taken as the s at the peak of figure left, in 
i bs- The extrapolated value, for < Q b s — > oo, is s c = 0.036 ± 0.001). The arrow indicates 
the theoretical value from Equation l(9]l (c) Susceptibility at s = s c = 0.036 635 37 
constant versus t Q b s - The critical divergence of the \k * s logarithmic. 



Let us now consider the system in the presence of an external magnetic field h, 
which couples to the spins of the chain according to Equation ©. We have studied 
the Ising Glauber chain with the magnetic field h close to the critical point. 

We remark that the magnetic susceptibility of the 2d Ising model is given by 

^- !S ^-Wit<^r l) . (1 3) 

1,3=1 

where G(ri,r 3 ) = ((Jicrj) — is the two point correlation function, so that Xj 

is static, and is a measure of the spatial correlations. In the critical isotherm, the 
magnetic susceptibility shows a critical divergence given by 

X T oc h-W> (14) 

being the critical exponent 5 = 15, thus 1 — 1/5 = 14/15 w 0.933. 

Following the definition of the generalized dynamic susceptibilities, we define 
the time-spatial magnetic susceptibilities as 



(M 2 ) - (M) 



2 1 P^ohs P^ohs 



kTNt ohs kTNt ohs 
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where now G{r i: rj, t', t") = (ai(t')(jj{t")) s — {<Ti(t')) s {<Tj(t")) s is the two point and two 
time correlation function. This susceptibility (% m ) has both spatial and temporal 
correlations. We can also define in the s-ensemble the static (or instantaneous) 
susceptibility in the same way as in Equation ((13|) , which is given by 

W = (™ 2 )s ~ HI =J_ST G (r- vt = 0) (16) 
X T kTN kTN ^ 1 n ' J ' T j ' { } 

where we have introduced r = \t" — t'\. 

We can apply the schema used for standard critical phenomena given by the 
dynamic scaling hypothesis, and obtain a relationship between the critical exponents 
corresponding to these susceptibilities 1123 II . So, in the thermodynamics limit, 
assuming a translational and temporal invariant correlation function G, and the 
existence of a diverging correlation length £, we obtain 

R 

and 

where we have introduced i? = | — rj \ . Here z is the dynamical critical exponent. 
For our system we have taken T = T c , that is the critical isotherm, at constant 
s = s c . Therefore the length correlation diverges according to the magnetic field. 
By assuming the usual Ising behaviour [24H, we have 

£oc|/i|- e ,with = z/(/3 + 7 )" 1 . (19) 



From the Equations (fT7|) and ((T9l> we recover the static result a h ^ 1 ^ <5 ' ) / whilst 
from the Equations (jl8|) and ((19|) , we find the dynamical critical behaviour given by 

Xm cx/i- A ,with A = ( 1 + vz)(/35)-\ (20) 

Figure H] shows the results of the Monte Carlo simulation of the Ising chain with 
a magnetic field for 5 x 10 6 generated trajectories (N = 64 and t obs = 300). A clear 
behaviour emerges for external fields of moderate intensity, that is 0.004 < h. We 
have found that the results for different times t b s = 160, 200, 300 (and different 
number of spins iV = 64, 100) are the same within error bars. For these values of 
the magnetic field there are no finite size effects. The stationary thermodynamic 
behaviour is observed. In other words, for weak field h < 0.004, there are finite 
size effects observed, caused by limited t obs and lattice size N. Any attempt to study 
the phase space region with a small external field (h < 0.004) would require extensive 
simulation efforts which are beyond the scope of this paper. 

Figure Hta) shows the static susceptibility x^)- We can see that a power law 
fits in the region of the magnetic field, 0.004 < h < 0.2, with a fitted exponent 
0.90 ± 0.05. For higher field (0.2 < h) there is a deviation of the power law regime, 
as expected far from the critical point. This value is consistent with the value of 
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the 2d Ising model given by the Equation (|T4~1) . On the other hand, figure Hfb) 
shows the susceptibility x m - The behaviour of the static susceptibility is qualitatively 
analogous to x m - However the power law fit within the same region of the magnetic 
field, 0.004 < h < 0.2, gives a different exponent, that is A = 1.50 ± 0.05. Using 
the relationship given by Equation (|20|) and assuming the Ising universality, we 
can conclude that our simulation gives us an estimation of the dynamical critical 
exponent, which is z = 1.06 ± 0.09. This value is compatible with z — 1. 

This results shows that through the application of the magnetic field to the Ising 
chain using the s-ensemble, the spatial correlations involved are similar to those 
present in the 2d Ising model with field near the critical point, while the response 
of the integrated magnetization is governed by a (greater) exponent dependent on 
the dynamical exponent z. 

In the system presented here, the static criticality of the 2d Ising model 
universality class appears in different ways. For the energy and activity criticality 
arises in the space-time plane (see figures |2] and |3|) and is uncovered by means of 
finite size effects, whereas for simulations with external magnetic field it appears 
only in the spatial axis (see figure lUja)) and in the thermodynamic limit. 

3.3. Dynamic first order transition 

We argued in section 13. 1[ that in the presence of an external field, the spontaneous 
broken symmetry should give rise to a first order phase transition in the 
ferromagnetic phase. So that, for isothermal conditions, a first order phase transition 
is expected for the biasing field s greater than the critical value s c (T ) (see figureUJ^c)). 
Illustrating the dynamic behaviour, including an external magnetic field analytically 
is non-trivial. This is due to the fact that the mapping of the master equation onto the 
quantum problem will give rise to highly non-linear terms. In the following we will 
address this problem numerically. 

We consider the parameter space {T, s, h} with focus on an isothermal plane 
(T = 1). In the h = plane, two phases are present; a paramagnetic phase for 
s < s c ps 0.037, and for s > s c ps 0.037 a ferromagnetic phase, as depicted in the 
schematic of figure [TJ Figure [5] shows the results of the Monte Carlos simulations 
including an external magnetic field, that is the magnetization in terms of h for 
different values of the biasing field s. In this case we have used a system of 
iV = 200 spins, and we generated 10 6 trajectories. In figure |5fa) a smooth curve 
m(h), in the vicinity of h = 0, for a low biasing field s = 0.02 can be seen, this is 
the paramagnetic phase. Whereas figure |5jc) shows a sharp transition, at h = 0, 
for a supercritical value of biasing the field s = 0.05, this is the ferromagnetic 
phase. For completeness figure |5{b) illustrates the behaviour at the critical point 
s c = 0.036 635 37. The behaviour of s is equivalent to that of the inverse temperature 
for the thermodynamic fluctuations of the Ising ferromagnet. Therefore, we expect 
that the field-driven first order phase transition ends at the critical point s = s c 
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Figure 4. (a) Susceptibility of the instantaneous magnetization corresponding to 
spatial fluctuations. The continuous line is a fit of a power law that gives exponent 
0.90 ± 0.05, in good agreement with the exponent of the magnetic susceptibility of 
the 2d Ising model 1 — 1/8 « 0.933 (b) Susceptibility of the integrated magnetization 
corresponding to spatio-temporal fluctuations. The continuous line is a fit of a power 
law that gives exponent 1.50 ± 0.05; 

(T — 1, h — 0) which is the limit of the ferromagnetic phase (s > s c ). 

In order to capture this behaviour further, hysteretic loops were studied. 
The simulations were started from a representative trajectory of the stationary 
distribution, then the external field is perturbed by h ± Ah (with Ah ~ 0.01). In 
order to reach a stationary state for the new value of h, n re \ ax steps in the TPS 
algorithm are taken for equilibration. Then measurements of the observables are 
taken for n ohs steps. Figure [6] shows the resulting loops of the magnetization versus 
the magnetic field. We used a system size N = 200 and n relax = 10 2 TPS moves 
were attempted, in order to achieve relaxation and the hysteresis was studied using 
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Figure 5. Stationary results for the magnetization of the Ising chain with magnetic 
field using the s-ensemble for different values of the constant biasing field s. The 
temperature T is constant (T = 1). (a) magnetization for the subcritical behaviour, 
that is in the paramagnetic phase s = 0.02 < s c , (b) results for the critical values of 
s = s c = 0.036 635 37, and (c) results obtained for the supercritical behaviour, that is 
in the ferromagnetic phase s = 0.05 > s c . 

n bs = 10 5 realizations. 

The top panel of figure |6] depicts the hysteresis loop with increasing supercritical 
s with a trajectory length of t bs = 60, that is s = 0.05 for (a), s = 0.07 for (b), and 
s = 0.08 for (c). The bottom panel of figured again shows the hysteresis loops, but 
in this case the observational time of the trajectory was increased to t G bs = 80. The 
values for s are retained (s = 0.05 for (d), s = 0.07 for (e), and s = 0.08 for (f)). Clearly, 
the area of the loops increases with t Q b s as well as with s. This is because as s increases 
the metastability of the system increases, therefore the loop area also increases. The 
behaviour with respect to t Q b s can be treated in terms of an increase in volume, which 
means the relaxation is slower and therefore an increased loop size is observed. This 
behaviour supports that the dashed line in the figure [He), and the surface below of 
the continuous line in the figure []Ja), are made of first order phase transition points. 

4. Kawasaki dynamics 

So far we have looked at dynamical models with non-conserved magnetization 
[fl4l H7| , which in the case of the classic thermodynamics, would serve as an order 
parameter. These are A type models or non-conserved order parameter [25 j. In the 
particular case of the Glauber dynamics it has been shown, and we have verified that 
in the dynamic exploration of this model, a line of critical points in the phase space 
can be found. We expect this kind of behaviour for model A type dynamics. 

On other hand, it is interesting to ask about the role of the the underlying 
dynamics on the observed dynamical phase diagram. In order to illuminate this point, 
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Figure 6. Magnetization versus the magnetic field in the form of hysteresis loops at 
constant temperature T = 1 around the first order transition at h = 0. Also the biasing 
field s is fixed (all of them supercritical values): for (a) and (d) s = 0.05, for (b) and 
(e) s = 0.07, and for (c) and (f) s = 0.08. In the top (bottom) panel we have used an 
observational time i b s = 60 (t b s = 80). The area of the loops increases with i b s , and 
with s. More details in the text. 

let us now consider a model belonging to the second class of dynamics, where the 
magnetization is conserved (conserved order parameter). We have realized Monte 
Carlo simulations using Kawaksaki dynamics. This dynamics conserves the overall 
magnetization and consists of exchanging neighbouring spins, rather than flipping 
them. The acceptance function for undergoing a spin exchange is chosen to be the 
Metropolis criterion 

P accept = Min(l,e- A ^ T ). (21) 

where AE is the increment in the energy due to the exchange. 

In figure [7| we show the results of the Monte Carlo simulations using the s- 
ensemble with Kawaksaki dynamics for the Ising chain. These simulation were 
realized at constant temperature T = 3 for a spin chain of size N = 64 and 
zero magnetization, i.e. equal populations of spins up and down. Figure ITfa) 
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Figure 7. Results of the Monte Carlo simulations using the s-ensemble with Kawasaki 
dynamics taking the temperature constant (T = 3). When the observational time i bs 
is varied, as indicated in (a), finite size effects are observed, (a) Activity versus the 
biasing field s. (b) susceptibility of the activity as a function of s. 



and figure IZtJb) show the activity and the susceptibility of the activity given 
by Equation (0, respectively versus the biasing field s for different t obs - For low 
values of s we can observe an active phase which becomes inactive for higher values 
of s. Also, when the observational time is increased, the jump becomes more abrupt 
and consequently its susceptibility shows an increasingly sharp peak. Again this is 
the sign of a phase transition and suggest a good finite size scaling behaviour through 
increasing the length scale, i.e. t ohs , in the same way was done in section|3l 





Figure 8. Results of the Monte Carlo simulations using the s-ensemble with Kawasaki 
dynamics taking the temperature constant (T = 3). When the observational time t a h s 
is varied, as indicated in (a), finite size effects are observed, (a) Energy versus the 
biasing field s. (b) susceptibility of the energy as a function of s. 



Figure [8] shows the behaviour of the integrated energy with respect to s and its 
susceptibility defined by Equation (fl2l) . This behaviour is qualitatively the same as 
the behaviour of the activity We can understand this similarity qualitatively The 
energy is proportional to number of domain walls (that is the number of the pair 
of first neighbours with different spins). Using the Kawasaki dynamics we require 
two neighbours sites with different spins, in order to make a spin flip, thus the same 
factor governs both the energy and activity. 
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The figures [7] and |8] suggest, that the observed transitions are first order 
transitions. Generally finite size scaling of a first order transition ||2T| predicts scaling 
in terms of the volume V = L D . According to this theory, for a finite volume V, we 
can define an effective critical point s c as value of s at which the susceptibilities have 
the maximum ■y/ max . In this way, the following scaling behaviour is expected 

Sc(tobs) = + AV~ l , (22) 

where is the thermodynamic limit s c (t ohs — > oo), and 

X max oc V . (23) 

In our case, the volume in the space-time is given by V = Nt y, s . By fixing the number 
of spins N, we can write out V = t ohs in Equations (|22|) and (f23)) . 

Figure |9] illustrates that the finite size analysis holds for the activity (the same 
behaviour is also observed for the energy). In figure Ufa) the location s c of the 
maximum x~k ^ s plotted for two Ising chains of length iV = 64 and iV = 100 against 
t~ b g. The effective critical value s c decreases with increasing t ohs and with increasing 
N. Using Equation (|22|) to calculate the extrapolated value to t obs — > oo we obtain for 
iV = 64 the value s c = 0.003 ± 0.001, and s c = 0.001 ± 0.001, for N = 100 (the same 
results are obtained from the energy analysis). 

Figure Htb) shows the height of the peaks drawn from the susceptibility of the 
activity versus the observational time. The continuous line in this figure is a power 



law given by the Equation (|23|). For the smaller system size the power law is not 
observed at all. For the larger size we can see that there is a narrow region, for 
intermediate values of t b s/ where the power law holds. These result shows that for 
small t bs there are boundary effects in the time of the trajectories, and for longer 
t obs there are boundary effects in the space of the spin chain. These rounding effects, 



make it difficult to observe the scaling behaviour given by Equation (|23|) for the sizes 
and observational times used in this paper. 

The previous analysis of figures Ufa) and Htb) indicates a first order phase 
transition, between an active phase (for s negative) and an inactive phase (for s 
positive) and a transition point at s c = 0. This space-time phase structure, of 
coexistence of active and inactive dynamical phases at s = 0, is similar to that of 
idealised kinetically constrained models 0.H]]. 

It is worth mentioning that the activity and the energy, in the inactive phase 
(see figures [Zta) and[8{a)) are ~ l/N due to finite size effects. In fact the dynamics 
implies that there always are, at least two domain walls, which separate the up and 
down pointing spin domains in order to conserve zero magnetization. Because of 
this constraint the minimum of the energy per site is given by e — e = 2/N, where 
eo is the unconstrained energy corresponding to the ground state. For increasing iV 
and increasing t ohs the energy in this phase tends to this characteristic value, and 
consequently the activity is of the same order. This behaviour shows that the inactive 
phase is a ferromagnet, where the lattice phase separates into two domains, one 
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Figure 9. Results of the Monte Carlo simulations using the s-ensemble with Kawasaki 
dynamics taking the temperature constant (T = 3). Black Squares are results obtained 
from the data of figure [7] for N = 64. For comparison, red triangles are data obtained 
from simulations with N = 100 analogous to figure [7] (data not showed for the sake 
of clarity), (a) The effective critical biasing field s c , defined as the location of the peak 
in the susceptibility x^, as a function of l/i bs- The straight lines are fits used for the 
extrapolations of s c to the limit t ^ s — > oo, by means of the scaling valid for first order 
transitions given by Equation The result is s c = 0.003 ±0.001 (s c = 0.001 ±0.001) 
for N = 64 (N = 100). (b) height of the peak of the susceptibilities of the activity for 
the same cases that in (a). The straight line is the expected power law cc t Q b s from the 
finite size theory for first order transitions. More details in the text. 
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Figure 10. Phase diagram of the Ising chain with Kawaksaki dynamics. The dashed 
line for s = represents points of a first order phase transition between the indicated 
phases. 



consisting of up spins and the other of down spins. Such state will have zero activity 
in the thermodynamic limit. 

We have performed the same analysis for other temperatures, and we conclude 
that this first order dynamical transition between a paramagnetic /active phase and 
a ferromagnetic /inactive phase at s c = is present for all temperatures considered 
(T = 1, 2, 3). At low temperatures, the accessibility of the transition becomes more 
difficult because the value of activity will be much smaller at lower temperatures and 
it will be zero in the inactive phase, so that the difference between the two phases 
is too small to resolve. In the limit T — > the discontinuity disappears and the 
transition ends in the point s = 0, T = 0. Once again, this is analogous to what 
occurs in kinetically constrained models ©H]]. 

Figure [10] presents the phase diagram that summarizes the result discussed 
above. This phase diagram for the Ising chain using Kawasaki dynamics should 
be compared to figured] illustrating the Glauber dynamics. These diagrams suggest a 
very different phase behaviour for Models A and B. Models B only exhibit a dynamic 
phase transition at s = between two dynamic asymmetric phases, one is an active 
paramagnet and the other is an inactive ferromagnet. This first order transition is 
driven by the biasing field s. For Models A, there is a surface of first order lines, 
between two symmetric ferromagnetic phases and it is driven by the magnetic field. 
The surface ends in a critical line, in the 2d Ising universality at s ^ 0. 



Thermodynamics of trajectories of Id Ising model 



22 



5. Conclusions 

In this paper we have considered how dynamical phase behaviour relates to 
thermodynamics in the context of the one-dimensional Ising model. We have studied 
s-ensembles of trajectories [2H3}|I1[9| by means of path sampling simulation methods 
for both Glauber and Kawasaki dynamics. We have both confirmed the analytic 



predictions of reference [17] for Glauber dynamics, and extended the space-time 
phase diagram to the case of non-zero magnetic field. Our results show an intimate 
relation between thermodynamic and dynamic phases, but this relation depends 
on the details of the dynamical rules. For both Glauber and Kawasaki dynamics 
we show that the paramagnet (which is the stable thermodynamic phase for all 
T > in one dimension) has a counterpart dynamical phase of high activity. 
While thermodynamically unstable, the static ferromagnetic phase corresponds 
dynamically to a phase of low activity, and by biasing ensembles of trajectories via 
the counting field s it is possible to go through a dynamic/non-equilibrium phase 
transitions between the active paramagnetic, and the low-activity ferromagnetic 
phases. In the case of Glauber dynamics, as predicted in |fT7|, this transition is second- 
order and in the 2d Ising model universality class, a fact that we verified numerically 
by detailed finite size (and time) scaling. Furthermore, with the inclusion of a 
magnetic field h we were able to show that the area below the line of critical points 
in the T — s plane at h = is actually a surface of first order transitions between low- 
activity ferromagnetic phases of positive or negative overall magnetization. This is a 
very interesting example of dimensional reduction: ensembles of trajectories of the Id 
Glauber Ising model have statistical properties similar to ensembles of configurations 
of the equilibrium 2d Ising model. 

We also considered the Id Ising model under Kawasaki dynamics. In this case 
the transition line between the active /paramagnet and the inactive/ ferromagnet lies 
on the s = axis, a situation analogous to that of kinetically constrained models 
of glasses B2HU. This can be understood by the fact that under Kawasaki dynamics 
spin-exchange leads to transitions in terms of domain walls where there is at least 
one domain wall present either before or after the transition: the presence of this 
persistent domain wall plays the role of a kinetic constraint. 

The aim of this paper was to shed light on the connection between structural 
or thermodynamic phase structure, and space-time or dynamic phase structure 
as revealed by the large-deviation method of trajectories. The natural next, and 
more challenging, step in this programme will be to perform a similar analysis 
on prototypical systems which display thermodynamic transitions at non-zero 
temperature, such as the 2d Ising model. 
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